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Abstract 

The general expectation that, in principle, the time-dependent density functional theory 
(TDDFT) be an exact formulation of the time-evolution of an interacting A^-electron system is 
critically reexamined. It is demonstrated that the previous TDDFT foundation, resting on four 
theorems by Runge and Gross (RG) [Phys. Rev. Lett. 52, 997(1984)], is invalid because undefined 
phase factors corrupt the RG action integral functional. Our finding confirms much of a previous 
analysis by van Leeuwen [Int. J. Mod. Phys. B 15, 1969(2001)]. 

To analyze the RG theorems and other aspects of TDDFT, an utmost simplification of the 
KS concept has been introduced, in which the ground-state density is obtained from a single KS 
equation for one spatial (spin-less) particle. The time-dependent (td) form of this radical Kohn- 
Sham (rKS) scheme, which has the same validity status as the ordinary KS version, has proved 
to be a valuable tool for analysis. The rKS concept is used to clarify also the alternative non- 
variational formulation of td KS theory. We argue that it is just a formal theory, allowing one to 
reproduce, but not predict the time-development of the exact density of the interacting A'^-electron 
system. 

Besides the issue of the formal exactness of TDDFT, it is shown that both the static and time- 
dependent KS linear response equations neglect the particle-particle (p-p) and hole-hole (h-h) 
matrix elements of the perturbing operator. For a local (multiplicative) operator this docs not lead 
to a loss of information due to a remarkable general property of local operators. Accordingly, no 
logical inconsistency arises with respect to DFT, because DFT requires any external potential to 
be local. For a general non-local operator the error resulting from the neglected matrix elements 
is of second order in the electronic repulsion. 
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I. INTRODUCTION 



Over the last decade, time-dependent density functional theory (TDDFT) has become an 
extremely popular method to compute electronic excitation energies and response properties 
of ever bigger molecules and clusters (see, for example, Refs. . The method and its 



foundations were already worked out in the nineteen eighties, primarily in papers by Zangwil 



and Soven , Runge and Gross (5|] , and Gross and Ko^ 



computer codes have been developed 



12 



m p]. More recently, various efficient 



and made available as parts of major 



quantum chemistry program packages. At present, one witnesses intense activities worldwide 
aiming both at the further development of methodological aspects and the computational 
efficiency of the codes. 

Besides the relatively modest computational expense, a major boost for the advancement 
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of the method has been the assurance that TDDFT is a formally exact theory 
that is, the TDDFT results would become exact if the exact time-dependent exchange- 
correlation (xc) potentials were available. In practice, of course, one always has to use 
approximate xc potentials, and therefore one has to be prepared for smaller or larger errors in 
the computational results. There is a widely held confidence that any problems encountered 
with the TDDFT method are only caused by imperfections of the underlying xc potentials, 
a belief prevailing even as some more severe problems became apparent, such as in the 



description of Rydberg excitations |16|, the treatment of extended vr systems [17|, the 



absence of double (and higher) excitatioris 
transfer (CT) excitation energies 



1^ liol . and the 1/R dependence of charge- 
|23|. These failures have triggered efforts to 



modify the xc potentials accordingly and thereby remedy the respective problems. Most of 
this work has been confined to the so-called adiabatic approximation, in which the time- 
dependence enters the xc functionals only via the time-dependent (td) density functions. 
But also the development of time- or energy- dependent xc functionals beyond the adiabatic 
approximation [l^ Q has been envisaged. 

On the other hand, TDDFT has never obtained a similarly accepted status of uncontested 
validity as the original (time-independent) density functional theory (DFT) developed by 
Hohenberg and Kohn (HK) 26 1 and by Kohn and Sham (KS) [27[. The foundations of 
TDDFT, as formulated by Runge and Gross (RG) in Ref. P], have been constructed largely 
in terms analogous to the HK and KS concepts of DFT. However, elusive notions such as 
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td v-representability and non-interacting v-representabilty were clearly in need of further 



mathematical clarification 
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More recently, the RG foundations of TDDFT were 



challenged by Rajagopal j3l|, van Leeuwen j32], and Harbola and Banerjee [33[, after it was 



realized that the kernel of the xc functional in the RG formulation violates causality 
A critical review of the RG action integral functionals by van Leeuwen j32, 36| revealed basic 
deficiencies. Presently, an alternative formulation of TDDFT, being essentially a KS-type 



approach without imp 



foundation 
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nng a variational principle of the HK type, is viewed as a valid 



In this paper we will take a new look at the foundations of TDDFT. For our review we 
use a simple analytical device, referred to as radical Kohn-Sham (rKS) formulation, which 
is as legitimate as the usual A^-electron KS theory. Not obscured by intricacies such as 
td v-representability etc., the rKS concept allows us to analyze both the static DFT and 
TDDFT in an utmost transparent way. What we find, confirms van Leeuwen's criticism of 
the RG foundation of TDDFT, but also proves the non- variational form of TDDFT illusory. 

An outline of the paper is as follows. The starting point of our study (Sec. 2) is the 
observation that the TDDFT equations, more specifically, the time-dependent Kohn-Sham 
(KS) linear response (LR) equations, neglect matrix elements of the perturbing (external) 
potential of the h-h or p-p type, where h and p refer to occupied (hole) and unocccupied 
(particle) KS orbitals, respectively. Because the exact linear response depends on all matrix 
elements, the TDDFT results appear to be deficient irrespective of the choice of the xc 
potential. The same situation arises in the case of a static (time-independent) perturbation, 
as is analyzed in Sec. 3. Here the problem would even challenge the well-founded (time- 
independent) DFT. The resolution of that puzzle in Sec. 3 is a very instructive confirmation 
of the logical consistency of DFT. In Sec. 4 we introduce the rKS concept, in which the GS 
density is not determined from the density of non-interacting electrons, but from a single 
KS equation for one (spin-less) particle. In Sec. 5 the rKS formulation is used to analyze 
the RG theorems and other aspects of TDDFT. A reader primarely interested in the issue 
of the validity of TDDFT might skip Sees. 2 and 3 and leap directly to Sec. 4. A summary 
of our results and some conclusions are given in the final Sec. 6. 
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II. COMPARISON OF EXACT AND KOHN-SHAM LINEAR RESPONSE 
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The TDDFT formalism has been presented in various ways in previous work 



where the reader is referred to for an overview and further details. Because TDDFT 



(in linear response form) is similar to the time-dependent Hartree-Fock (TDHF) or random- 
phase approximation (RPA) 



4G 
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42| . it is rewarding to consult also previous TDHF 
derivations (see, for example, Ring and Schuck j^^). A few basic notions pertinent to the 
ensuing discussion will be given in the following. 

For an A^-electron system (atom or molecule) having a non-degenerate ground state |^o) 
the (exact) ground state (GS) density matrix, 7, is given by 

Ipg = (^olcjcpl^o) (1) 

Here the second- quantized operators c],(cp) are associated with one-particle states (spin- 
orbitals) (pp. As a particular choice, we will consider the KS orbitals arising from the KS 
one-particle equations associated with the ground state of the system under consideration, 

/i^^0,(r, s) = {-iV' + vir) + J[p]{r) + t;,e[p](r)}0,(r, s) = e^r, s) (2) 

Here r and s denote spatial and spin variables, respectively, v{r) is the one-particle operator 
for the electron-nuclei interaction, J[p](r) is the Coulomb operator, and Vxc[p]{f) is the KS 
exchange- correlation potential. The exact GS density function, 7(r), is obtained from the 
density matrix elements according to 

p,q s 

By contrast to the exact density matrix, the KS density matrix, p, is derived from the KS 
determinant 

1$^^) = 101... 0^1 (4) 

according to 

Pp. = {^r\clc,\^r) (5) 

In the KS orbital representation assumed here the KS density matrix assumes the simple 
diagonal form 

Ppq = SpqUp (6) 
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where = 0, 1 denote KS occupation numbers. The KS density function, 

p,q s k s 

is devised to reproduce the exact density function, that is, 7(r) = p(r), provided the correct 
exchange- correlation potential is used in Eq. Whereas, at least in principle, the exact 
and KS density functions are identical, the density matrices necessarily must differ. As is 
well recognized (see, for example, 3, ^^); the two entities differ with respect to a basic 
property: the KS density matrix, deriving from a single determinantal wave function, is 
idempotent, that is, = p, whereas the exact density matrix is not, 7^ 7^ 7. 

To discuss the linear response (LR) theory let us consider an additional time- dependent 
(td) external ( "driving" ) potential of the form 

u = df{t) (8) 

where d = d{r) is a local (multiplicative) operator and /(t) is a scalar time dependent 
function (vanishing for t < 0). In the notation of second quantization, the corresponding 
iV-electron operator, D = d{i), can be written as 

D = Y^ drsclcs (9) 

r,s 

where drs = (0rMl0s) denote the one-particle matrix elements of d. Now let us distinguish 
particle-hole {p-h) and h-p matrix elements, dak-, dka, from p-p and h-h elements, dab, dki- Here 
and in the following we use the notation in which the subscripts a,b,c, . . . and k, . . . 
denote unoccupied (virtual) and occupied KS orbitals, respectively, while the subscripts 
p,q,r, . . . refer to the general case. As will be discussed below, the KS linear response 
contribution to the density depends only on the p-h (and h-p) matrix elements of the driving 
potential, whereas the exact linear response contribution is a linear function of all matrix 
elements. 

Let us first inspect the exact case. Upon Fourier transformation the linear response of 
the exact density matrix can be written as (see, for example, {4^ ) 

Here |\E'„) and denote excited energy eigenstates and eigenvalues of the original (undis- 
turbed) Hamiltonian H; the complex infinitesimal irj is required for the definiteness of the 



Fourier transforms between the time and energy domain. For the special operator D con- 
sidered here, the transition moments appearing in the numerators on the right-hand side of 
Eq. ()10|) take on the form 

(*„|D|^0) = Y,dp<i{'^n\clc,\^o) (11) 

and it is obvious that the p-p and h-h contributions, (^„|c]jCb|\I/o) and (\I/„|c|.q|\E'o)) need 
not vanish. Using many-body perturbation theory (MBPT) for |\E'o) and \^n) based on the 
famihar M0ller-Plesset decomposition of H (and adopting for a moment Hartree-Fock (HF) 
one-particle states) one may readily establish that non-vanishing contributions appear for 
the first time in second order. For example, one finds 

(^„|clc,|^o) = 0(2) (12) 

for single excitations deriving from the HF configurations c\cj |$^^). An analogous result 
is found in the case of the h-h amplitudes. Moreover, the exact response comprises contribu- 
tions arising from double (and higher) excitations, the corresponding p-p and h-h amplitudes 
being here even of first order. 

Now let us turn to the KS response theory. In the formulation given by Gross and 
Kohn [6] (adopting here a slightly deviating notation) the linear response to the KS density 
function is given by 

5,(rM-fxiry;.H"(r',.W (13) 

Here 

I < \ M'r.s)<j)q{r,s)(P*{r',s')(t)p{r',s') 
X(r, r-u) = Y. - .-e, + e, + ^6 ^^'^ 

pq s,s' ^ ^ 

is referred to as the KS density-density response function and 

vl^^{r, u) = u{r, uj) + J[Sp\ (r) + Sv^^r, u) (15) 

is the first-order effective potential comprising the (Fourier transformed) external perturba- 
tion of Eq. (jH)), the (first-order) change of the Coulomb potential, J[Sp]{r), and of the xc 
potential, 6vxc{f,uj), the latter two contributions being linear expressions in 6p{r,uj). 

Inserting the first part of v\^^{r,uj), that is, the "driving" potential, u{r,uj), in the rhs 
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of Eq. (O yields 

X{r,r;u;)u{r ,Lj)dr =yy{ ■ dak (16) 

J ^ \uj - ea + ek + id 



a,k s 

(f)l{r,s)Mr,s) 
if + ea - ek + iS 



dka /(t^) 



Obviously, here the p-p and h-h matrix elements of d have been projected out and only p-h 
(and h-p) matrix elements, dak, enter (as inhomogeneities) the linear KS response equations. 
As we have seen, the exact linear response to the density function, 6'~f{r,u), which may be 
written in a form analogous to Eq. (fT!?jl . 

S-,ir,.) = jxiry:.Hr'..)ir' (17) 



is a linear function of all matrix e 
density-density response function 



^ments, dpg, of the perturbing potential. Here the exact 
x(r, r'; u), is related to ^7 of Eq. (|Tn|) according to 



X{r, r'- u)d{r')dr' = J2Y1 ^^mi^'W^^ ')Mr, s) (18) 

p,q s 

This result evokes the question if equating 6p{r, uj) and 5'-f{r, uj) is permitted at all. It seems 
that in the KS linear response equations the information associated with the p-p and h-h 
matrix elements of the external (driving) potential is lost and, thus, these equations have to 
be viewed as an approximation even in the case of an exact exchange correlation potential. 
Let us inspect the situation in the more general and transparent matrix formulation of 



TDDFT (see, for example. 



12, ll^). Here the KS response equations are written in the form 



of a matrix commutator relation, 

ij6p= [h, 6p] + [6h, p] + [d, p] (19) 

where h and d denote the matrix representations of the (unperturbed) KS Hamiltonian and 
the perturbing potential (time-independent part), respectively, and 6h is the change of the 
KS Hamiltonian linear in 6p. Note that due to the form of p the commutator [d, p] on the 
rhs of Eq. (fT^ projects out the p-p and h-h matrix elements of d. Arranging the p-h and 
h-p matrix elements of Sp and of d in columns (vectors), 

*P=| ^-"V I (20) 



the linear response equation for 5p takes on the famihar RPA form, 

5p = d (21) 




The elements of the matrices A and -B, being related to the functional derivatives of the 
Coulomb and xc potentials of the KS Hamiltonian, have been specified elsewhere (see, for 
example, Ref. [3], and Sec. 3). The information on the perturbing external potential enters 
the response equations ^2\\ji only via the vector d. Thus, it is manifest that only the p-h 
and h-p matrix elements of the perturbing potential come into play. Does the neglect of the 
h-h and p-p matrix elements of the perturbing operator mean that TDDFT is not formally 
exact? The same problem occurs in the case of a time-independent (static) perturbation, 
and here it would even challenge the logical consistency of DFT itself, more specifically, the 
universality of the HK energy functionals. In the next Sec. 3 we will consider the static case 
and see how the apparent contradiction to the universality of the DFT functionals can be 
resolved. 



III. TIME-INDEPENDENT KOHN-SHAM RESPONSE THEORY 

In this section we consider the problem of the loss of the h-h and p-p matrix elements in 
the simpler static case of a (small) time-independent external perturbation, u. 

The linear KS response equations are obtained here as a special case [u = 0) of the more 
general time-dependent equations (fT^ : 

[h,6p] + [6h,p] + [u,p] = (22) 

Likewise, these equations can be deduced via first-order perturbation theory for the ground 
state KS orbital, also referred to as coupled-perturbed Kohn-Sham (CPKS) theory (see 
Casida^^ and references therein). As above h and u denote the matrix representations of 
the unperturbed KS Hamitonian and the perturbing potential, respectively. The KS density 
matrix, p, associated with the unperturbed ground state is diagonal, ppq = SpqUp] 6p denotes 
the first-order change in the KS density matrix. Finally, 6h is the matrix representation of 
the linear change of the KS Hamiltonian, 

6h = J[6p]{r) + 6v^,{r) (23) 
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A basic assumption of the KS linear response theory is that the xc part of 5h can be expanded 
according to 



SVr.r. 



v.Ap + M - ^Ap] = I ^-^J^^Pi'^'^dr' + OiSp') (24) 



in terms of 6p (and, possibly, gradients, V5p, and higher derivatives). Here 6p is related to 
the first-order density matrix, 6p, according to 

^Pir) = XI 5Z s)5pak + 0a(^> s)0fc(r, s)5pka} (25) 

a,k s 

Proceeding in the usual way, one now may evaluate the p-h and h-p matrix elements of 

6h, 

{5h),, = j 5h{r)ct>;{r)<P,{r)dr (26) 
which leads to linear expressions in the density matrix elements 5prs'- 

{5h)pq = X Mpg^rs Sprs (27) 
rs 

Here the index pairs, (pq) or (rs), are either of p-h or h-p type. Finally, introducing 6h in 
that form in Eq. ()22|1 one arrives at the desired KS linear response equations, reading in 
matrix form analogous to Eqs. 



A B 

B* A* 



6p = -u (28) 



Here, the matrix elements of A and B are given by 

Aak,bl = (Ca — ^k)SabSkl + Mak,bh Bak,lb = Mak,lb (29) 

As in Eq. ()20|). the p-h and h-p matrix elements of Sp and u are arranged to form column 
vectors, 

lp=(^'A, M-l^A (30) 

\ ^Php J \^hp j 

As in the td case, the perturbation enters the set of linear equations via the inhomogeneity 
vector, tt, in which the p-p and h-h matrix elements of u are absent. 

It appears that we are facing a paradox here: on the one hand, we have just applied 
a valid first-order perturbation theory to the density within the KS framework arriving 
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at a seemingly deficient result; on the other hand, the full solution of the KS eigenvalue 
problem for the perturbed Hamiltonian must yield the exact density, so that also the result 
of first-order perturbation theory cannot be incorrect. The answer to this puzzle is that no 
information on the perturbing potential is lost in the absent p-p and h-h matrix elements 
provided that the perturbation is a local (multiplicative) operator [45]. This is due to a 
remarkable, though apparently not widely known property of local operators which may be 
stated as follows: 

Theorem 1. A local operator, v = v{r), is uniquely determined up to a constant by its p-h 
(and h-p) matrix elements with respect to a complete one-particle basis and an arbitrary 
partitioning of that basis into occupied (hole) and unoccupied (particle) one-particle states. 
A simple proof of this assertion is given in the Appendix. An interesting aspect here is that 
the proof assures merely the uniqueness (up to a constant) of the local operator but does 
not offer a way to reproduce the operator from its p-h matrix elements. It seems that for 
such a reconstruction one needs one of the diagonal blocks, that is, either h-h or p-p, in 
addition to the p-h block. Thus, the logical status of theorem 1 resembles that of the HK 
and KS theorems which prove the existence of universal xc functionals but do not provide 
for means to construct the functionals. 

Theorem 1 assures that the loss of the p-p and h-h matrix elements in the CPKS equations 
is no contradiction to the formal exactness of the theory, provided that the external potentials 
are local (multiplicative). Clearly, this observation applies also to the td KS linear response 
considered in Sec. 2. The restriction to local external (one-particle) potentials is a basic 
and well understood consistency requirement of DFT. The universality of the kinetic energy 
and xc functionals hinges on the condition that the external potential functionals are of the 
form 

E,[p] = J v{r)p{r)dr (31) 

It should be recalled that quantum theory is essentially non-local, and many physically 
important interactions are not of the local type. For example, the interaction of electrons 
with an electromagnetic field involves the momentum operators, p. = —iVj. It is common 
practice, to apply the usual CPKS and TDDFT methods also for non-local external poten- 



tials (see, for example, 42, ^^^^ "^^^e one should be aware that the loss of the 

p-p and h-h matrix elements of the external operator introduces indeed an error beyond 
the approximation for the functional, which is of second order in the Coulomb repulsion. 
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In the td KS linear response this error affects only the transition moments. In particular, 
it destroys the equivalence between the so-called length and velocity forms of the transi- 
tion moments, because the former is associated with a local operator and the latter with a 
non-local operator. 

In principle, DFT can be extended to account for non-local external potentials as well. 
For this purpose the non-local potential must be incorporated "a priori" into the HK and 
KS formalism, that is, the HK and KS functionals have to be defined from the outset for the 
A^-electron system under the action of the external potential. This would lead to modified 
functionals being now specific to the considered non-local potential. In other words, the 
functionals would depend on the respective non-local external potentials. An important 
example for the necessity to deal with non-local operators is the presence of magnetic fields. 
As a systematic approach referred to as current density functional theory (CDFT) one here 
considers functionals that depend not only on the density but also on the current de nsity 
|4^0,0|. As another possibilty of dealing with non-local external potentials, Gilbert 
and Levy [0|] have considered density-matrix dependent functionals. 

Let us briefly inspect how the CPKS equations will change if the the xc potential depends 
directly on a non-local perturbing potential, u. Obviously, this would lead to an additional 
contribution to 6h of the form 

Sv:^c = w[p,u]{r) (32) 

and, thus, to another inhomogeneity term in the linear response equation p8|) . Here p is 
the unperturbed GS density. Because the additional inhomogeneity contribution depends 
on u, the full information on u can be restored, reconciling the (first-order) result of the KS 
linear response with the exact result. 

As a more general aspect, the non-local potential problem shows that the CPKS equa- 
tions, while justified as a valid first-order perturbation theory for the KS approach to deter- 
mine the ground state density of the system plus perturbation, may not be seen as physical 
response equations for the interacting A^-electron system (in that case they should apply 
also to non-local perturbations). This admonishes us to be wary of the prospect that the td 
KS equations can describe the time evolution of the system in response to a td perturbation. 
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IV. A RADICAL KOHN-SHAM VERSION 



The KS formulation is a clever way to transform the problem of finding the density 
minimizing the HK functional into the determination of the ground state of an associated 
non-interacting N-particle system. While providing a good basis for practical computational 
schemes, the usual KS formulation still does not achieve its full theoretical potential. In 
fact, one may proceed to a radical KS approach, in which the mapping of the exact density 
is not onto that of a non-interacting N-particle system but rather to the density of a single 
particle. Whereas such a radical KS formulation will be less suitable as a starting point 
for the approximate treatment of the exact ground state density, it may serve as a valuable 
analytical tool to clarify various aspects of DFT and, in particular, TDDFT. It should be 
noted that the idea of such an obvious extension of the usual KS approach is not new, 

n 

though apparently little known. Already in 1984 it was used by Levy et al. j5J| to discuss 
asymptotic properties of the xc potential. 

In the usual (N-particle) KS formulation the kinetic energy contribution, T[p], to the HK 
functional is substituted by the kinetic energy functional 

Ts[p] = J2[ ^^*(^' s)hlV'mr, s) dr (33) 

i,s 

of a non-interacting A^-particle system, the density being obtained according to 

k s 

as the density function associated with the Slater determinant |$) = \1IJ1ip2 ■ ■ - iPnI of or- 
thonormal orbitals = 1, . . . , N. As is well-known, the functional (|33|) can be made 

I (see also Parr and 



unambiguous using the Levy constrained search (LCS) definition 
Yang |55|): 

Ts[p]=mmm-lj2'^m (35) 

where $ ^ p indicates that the search is over all Slater determinants yielding the given 
density p. 

The deviation between the exact and the KS kinetic energy, T[p] — T5[p], is accounted 
for in the KS exchange-correlation functional, 

E,Jp] = T[p] - Ts\p\ + Kjp] - J\P\ (36) 
13 



so that the original HK energy functional can be written as 



Fhk[p\ = Ts[p\ + J\p\ + E,,[p\ (37) 

Here, Vee[p\ and J[p\ denote the full and classical electronic repulsion functional, respec- 
tively. Now the task of finding the minimum of the total energy functional, 

E[p] = FhkIp] + j p{r)v{r) dr = E[p{iji}] (38) 

under the constraint J p{r)dr = N can be performed in orbital space (see Parr and Yang 
[ssjl). The corresponding variational procedure yields the well-known KS equations for the 
ground-state of a system of non-interacting electrons moving in the external potential 

y^//[p](r) = vir) + J[p]{r) + v^,[p]ir) (39) 

where 

is the KS exchange-correlation potential. More precisely, v'^^f[p]{r) is a potential-functional, 
and Eqs. (|39l4m) have to be solved self-consistently. Self-consistency will be attained for 
the exact ground-state density, po, where the KS equations with the potential v'^^^lpolir) 
reproduce po- It should be noted that the consistency of the KS orbitals with the LCS 
requirement must be assured (see Levy and Perdew IsG^]). 

As this rigorous derivation of the KS equations shows, there is nothing that would compel 
a density representation associated with non-interacting electrons: any number of non- 
interacting electrons will be permissible, even A^ = 1. Indeed, we will demonstrate in the 
following how the entire line of arguments can readily be transferred to the representation 
of the density in terms of a single (spinless) particle. 

Obviously, any N-electron (ground-state) density function, p(r), can be represented by a 
one-particle wave function (orbital) according to 

p(r) = Ar|0(r)|2 (41) 

where 

<P{r) = {P^f' (42) 
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Such a representation is unique as long as p(r) > and p(r) 7^ for finite values of \r\. 
Obviously, this defines directly (that is, without invoking the concept of non-interacting v- 
representability) a 1-1 mapping of density functions and (real) orbitals. Next we can define 
a corresponding kinetic energy functional: 

fs[p] = N J 0(r)(-iV2)0(r) dr (43) 

Since the (real) orbital 0(r) is uniquely defined by p{r), so is the kinetic energy. 

It should be noted here that this definition is consistent with the Levy constrained search 
procedure. The general form of an orbital reproducing the density p{r) according to Eq. ()41|) 
reads 

^(r) = e'''^''U{r) (44) 
where k{r) is a real function. Clearly, the kinetic energy of ip, 

(^1 - iVV) = I / (VA:(r))^0(r)2 dr + (0| - |VV) (45) 

is larger than the kinetic energy of the real orbital 0, if k{r) 7^ const. This means that 
the orbital minimizing the kinetic energy functional for a given density is (up to a constant 
phase) a real function. As a consequence, Eq. (jl^ relates densities and orbitals, and the 
kinetic energy functional (j43|l is uniquely defined at the orbital level. We may elaborate that 
point somewhat further by considering a system where the KS orbital cannot be chosen real, 
e.g. in the presence of an external magnetic field. Clearly, an orbital of the general form (j44p 
is not determined by the density alone. In addition, one has to take into account the current 
density, j = 0^V/c, in order to obtain a unique definition of kinetic energy functional, now 
being a functional of both p and j, at the orbital level. That is why a current density version 
of DFT must be used in the case of magnetic fields. 

The next step is to introduce a correspondingly modified xc functional, 

E^M = T[P] - Ts[p] + Veeip] - J[p] (46) 

so that the functional for the total energy can be written as 

E[p] = fs[p] + J[p] + 4c[p] + J p{r)v{T) dr (47) 

As in the usual KS approach, the variational search for the minimum of E[p] under the 
constraint j p{r)dr = N can equivalently be effected by a search in the space of (normalized) 
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orbitals (f){r). The variation of £'[p{0}] with respect to (f){r) via 

p{r) = N\<P{r)\' 



(48) 



leads to the single KS equation 

{-iV + v'"[p]ir)}(l){r) = e0(r) (49) 
for the ground-state of a single particle moving in the effective potential 

^^//[p](r) = v{r) + J[p]{r) + i).c[p](r) (50) 
where Vxc[p]ir) is the modified xc potential deriving from Exc[p], 

Obviously, the single GS KS orbital has no direct physical meaning. It may be viewed as a 
kind of a mean orbital averaged over the N/2 spatial KS orbitals of the usual approach. 

As the usual KS approach, the radical Kohn-Sham (rKS) formulation, established by 
Eqs. ()48|49|) is, in principle, exact. That is, one would obtain the exact ground-state density 
of the interacting N-electron system provided the exact energy functional were available. 
Of course, the usual N-electron KS formulation will be a better starting point for the use 
of approximative functionals, simply because its expression for the kinetic energy, Eq. (jHSI), 
will give a better approximation to the full kinetic energy than the mean one-orbital term 
of Eq. (jlHj) . The actual benefit of the rKS variant is its potential as an analytical tool, and 
in the ensuing Sec. 5 we will use that tool to examine the foundations of TDDFT. 

Let us emphasize once again that Vxc[p\{i') is a potential- functional and Eqs. ()49l50p have 
to be solved self-consistently to yield the exact GS density, Po(^)- As a consequence of the 
simple structure of the rKS equation ()49|1 . the xc potential for the exact GS density, po{r), 
can be expressed according to 

v.c[po]{r) = ^ i— ^ V^yM^-t;(r) - J[po]ir) + e (52) 
Eq. (jK^ has been used to study features of the exact KS xc potential, such as the asymptotic 



behaviour 



3- 



A similar equation, arising in the ordinary 



K^S treatment of 2-electron 



systems, was used to characterize 2-electron KS xc functionals 



57 



58, 
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While the rKS version introduced above is the simplest possible KS-type approach, other 
variants are conceivable in which the non-interacting system consists of M = 2 or more 
electrons (M might even be larger than N.) For example, in the case M = 2 any (reasonable) 
density can be derived from the KS determinant 

l^o'''') = l0Oa0O/3| (53) 

for two non-interacting spin-| particles in the spin-orbitals, = (t^o{f)X'y{.s), 7 = 
where the spatial orbital is given by 

Mr) = (64) 

Whereas the rKS formulation is purely spatial, spin degrees-of-freedom come into play in 
these M-electron KS variants for M > 2. 



V. REVIEW OF TIME-DEPENDENT DFT 

A. Time dependent radical Kohn-Sham theory 

Having established the rKS formulation for the static case, we may now use this tool to 
analyze td density functional theory. 

Let us assume a td external potential, U(t) = '^u{ri,t), vanishing for t < 0, and let the 
system be in its (unperturbed) ground state at t = 0. The solution of the td A^-electron 
Schrodinger equation, 

gives rise to an associated exact td density function, p = p{r,t) with p(r',0) = po- As in 
the static case, the time development of the exact density can be assigned to a td orbital by 
generalizing Eq. 

0(r',t) = (^)V^ (56) 

This is trivial. The non-trivial issue is, of course, whether one can establish a Schrodinger- 
type equation at the single-orbital level that would allow one to predict the time-development 
of the exact density. Because inevitably any (non-stationary) wave function evolving accord- 
ing to a td Schrodinger equation picks up a time- and space-dependent phase, the orbital 
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must be written in the general form, 

^/;(r,t) = e*'=("'*V(r',t) (57) 

where k{r,t) is a real- valued phase function, and (p{r,t) is given by Eq. 

Supposing the RG theorems valid in their original form, they will apply as well to the 
rKS formulation. Then there is a single td KS equation of the form, 

^§-^Hr,t) = {-^W^ + vir) + J[p](r) + ^(r, t) + fi,.,[p](r, t)}V^(r, t) (58) 

governing the time evolution of ^jir^t) and, thus, of p{r,t) = A^|-?/'(r, t)p. Here 

is the td xc potential associated with the rKS modification of the RG td xc functional, 
^xc[p(^)]- The RG theorems assure that such a td xc potential exists, so that, in principle, 
the time development of the density can be determined exactly via Eq. ()58p. In practice, of 
course, one has to resort to approximations such as the widely used adiabatic local density 
approximation (ALDA). Here one uses the ordinary DFT xc potentials, 

a^c[p]ir,t) = v^c[pir,t)] (60) 

depending on time only via the time dependence of the density function, p = p{r,t). As 
above, the tilde indicates the rKS form of these quantities. 

In the rKS version the RG theorems suggest that one can, at least in principle, condense 
the full A^-electron td Schrodinger equation into a one-orbital td KS equation. Can this be 
true? As a step towards an answer let us inspect how the fourth RG theorem, establishing 
an analogy to the KS concept in the time-independent DFT, works in the rKS case. 



B. The Runge-Gross theorems 

The KS equations have been invented as a means for determining the minimum of the 
HK energy functional and thus the exact ground- state density of the interacting iV-electron 
system. In TDDFT the role of the KS equations is daringly expanded: their time-dependend 
form is believed to govern, at least in principle, also the exact time evolution of the density 
of the interacting N-electron system. The basis for that claim has been laid in a series of 

18 



four theorems in a famous article by Runge and Gross P|, in the following referred to as 
RG. Let us critically review their arguments. 

The first theorem (RGl) is the td analogue of the first HK theorem. It establishes a 
one-to-one correspondence between td density functions, p(r,t), and td external potentials, 
u[p]{r,t), which, in turn, via the td Schrodinger equation, 

d 

i- 



dt 



^[p]it)=(^H + U[p]it))^p]it) (61) 



determine the exact td N-electron wave functions, ^E'[p](t) (up to a purely time-dependent 
phase). 

The third theorem (RG3) is the analogue to the second HK theorem. Instead of the HK 
energy functional, one considers the action integral defined according to 

A[p] = J\t mp]m^^ - HMpm (62) 

We may leave any problems in this definition (see Ref. j3|) at that and go on further 
to the fourth theorem (RG4). In analogy to the ordinary KS approach, one introduces a 
kinetic-energy action functional, 

So[p] = jymp]m^^-fmp]{t)) (63) 

for non-interacting particles. Here it is supposed that for a given td density function, p(r, t), 
there exists a unique state (Slater determinant), <I>[p](t), of the non- interacting electron 
system. The functional 5*0 [p] is defined in analogy to the full kinetic energy action functional, 

s[p] = iymp]m^^-fmp]{t)) (64) 

for the original interacting electron system. As in the time-independent KS approach, So[p] 
replaces S[p], the remainder, S[p] — So[p], being transferred into the exchange-correlation 
part, Axc[p], of the full action functional (jHS))- Everything seems to be completely analogous 
to the time-independent case. 

However, there is a problem, clearly to be seen in the focus of the rKS formulation. 
Here the non-interacting state, $[p](t), becomes a one-particle state of the general form of 
Eq. (gH), 



19 



so that the So functional reads 



'0 

So[p] = dt I drr{r,m^^ + l^'Mr,t) (65) 

While the modulus of ip{r,t) is completely determined by the density p{r,t), the phase 
function k{r, t) is not. Clearly, the value of Sq[p\ depends manifestly on this phase function, 
but there is no way of determining it from the given density. This means that the functional 
5*0 is ill-defined at the orbital level. There are (infinitely) many orbitals for a given density, 
each giving a different value for the 5*0 functional. 

Let us consider the latter argument in somewhat greater detail. Inserting the form ()57|) 
of the orbital in the integrand of the Sq[p\ functional one readily obtains 

(V^l^l + iV^lV^) = (0||V2|0) - (0|(Vfc)2|0) - (0|A:|0) (66) 

This means that besides the density here also the gradient of the phase function, VA;(r,t), 
and the time derivative, k{r,t), is needed. Indeed, the latter information can be derived 
from the density to a certain, yet insufficient extent. Obviously, the orbital not only is to 
reproduce the density, but also to fulfill a td Schrodinger equation (SE) of the form 

z^^(r, t) = {-iV^ + w{r, t)}^(r, t) (67) 

where w{r, t) is a local td potential yet to be determined. Therefore, the continuity equation 

d 
Jt 

applies to the orbital, where the current density is given by 

j = (p'^Vk (69) 

As a consequence, it is possible to determine V/c from p and p, respectively. (A mathematical 
complication may arise here due to the requirement that V x (j/0^) must vanish.) Further, 
if Vfc is given (at any time), then also k{r,t) is determined, though only up to a pure 
time-dependent function, a{t). But as Eq. clearly shows, the latter indefiniteness of the 
phase function prevents the 5*0 functional to become well-defined. The time integral J a{t)dt 
on the r.h.s. of Eq. (|66p leads to a completely undetermined constant in the definition of 
the 5*0 functional. Note that this does not mean just a uniform shift of the 5*0 values, which, 
of course, would drop out in a variational search for stationary points. 
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j/ + yj-<> (68) 



So far we have not specified the local potential w{r,t) in Eq. (jHIj), but only assumed 
that such a potential exists, e.g., as a consequence of the first RG theorem (RGI) applied 
to the non- interacting KS system (of one electron). But in the rKS version, the RGI result 
can be obtained in a direct way (thereby proving the one-orbital td v-representability of 
any "reasonable" density), allowing us even to give an explicit expression for w{r,t). This 
is achieved by inserting the ansatz (j37j) in the SE ()67p. Separating the real and imaginary 
parts yields the following two equations: 

w{r,t) = ^-l{Vkr-k (70) 

and 

<j) + VkV(j) + lV^k(j) = (71) 

Obviously, the latter equation reproduces the continuity equation (jHHjl . whereas the former 
gives an explicit expression for the local td potential, w{r,t). Since both Vfc and k result 
from p as discussed above, w{r, t) is determined by p up to a purely time- dependent function, 
namely a{t). This shows that even the explicit form of the potential is of no avail to 
determine a{t). Even if the value of a(t), was given (or fixed) at an initial time, it cannot be 
determined for later times by solving the td SE ()(i7|l due to the corresponding indefiniteness 
of w{r, t). 

Thus, the rKS formulation inevitably points our view on the problem of the undetermined 
purely time-dependent phase functions corrupting the RG action integral functionals. The 
phase problem arises not only in the functionals of the non-interacting KS system but already 
in the functional (j62|l for the original interacting A^-electron system. When one consults the 
RG paper with regard to this issue, one finds that the phase problem, being discussed in 
the beginning of the paper, gets lost in the matrix element {^{t)\i^—T—W—V{t)\(^{t)) after 
Eq. (11). Here V{t) is the external td potential of the physical system under consideration. 
This potential does not contain a function C (t) that would cancel the time derivative of the 
phase function in $(t). Erroneously, RG argue here as if this potential was the potential 
V{t) according to the RGI theorem, that is, the potential invoked in the td SE for the 
A^-electron wave function, $(t), corresponding to the considered density. (Later TDDFT 
papers and virtually all review articles inconspicuously leap over the phase problem in the 
action integral functionals.) 

Whereas there is still wide-spread confidence in the RG foundations of TDDFT, their 
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breakdown due to the phase problem 



Leeuwen aheady several years ago 
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has been clearly analyzed and expressed by R. van 
In his 2001 review article js^ he draws the 
following conclusion: " We therefore conclude that time- dependent density-functional theory 
can not be based on the usual variational principle, and indeed attempts to do so have led 
to paradoxes." Let us note that besides the phase problem, van Leeuwen also analyzed 
correctly the non-stationarity of the RG action integral functionals, exposing another fault 
line in the original RG argumentation. But why have van Leeuwen's revelations not triggered 
stronger shock waves in the TDDFT community and beyond? Apparently because by the 
time of van Leeuwen's analysis the leading actors in the field had come to the conclusion 
that the KS equations could be established directly without the necessity of resorting to 
a variational principle jsj, 6^. Sharing that conviction, van Leeuwen communicated the 
reassuring message that TDDF is valid, though only in a new shape featuring the so-called 
Keldysh Green's function technique In the ensuing subsection 5.C, we will have a closer 
look at the non-variational KS theory. 



C. Kohn-Sham equations without a variational principle? 

In the derivation of the static KS equations three elements are essential: i) a universal 
energy functional (HKI); ii) a variational principle for the exact ground-state density (HKII); 
and Hi) a functional for the kinetic energy of non- interacting electrons defined at the orbital 
level (KS). Runge and Gross have pursued a strictly analogous approach in order to establish 
a basis for TDDFT. As first analyzed by van Leeuwen and corroborated here, this endeavor 
must be viewed as foundered in each of the three essentials. 

But is there a different route to establishing td KS equations? Within the TDDFT 
community, the generally accepted view is that this is the case. Indeed, the first RG theorem 
offers a shortcut to KS-type equations. Applying RGI to the case of non-interacting 
electrons, one can establish the mapping 

pir,t)-^w[p]ir,t) (72) 

so that the td KS-type equations 

t^Mr, t) = {-|V' + w[p]{r, t)}^,(r, t), j = l,...,N (73) 
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allow one to calculate the density p{r,t) from the orbitals ipj{r,t). The KS potential in 
Eq. ()73|) can be written in a more familiar form, 

w[p]{r,t) = Ve,t{r,t) + J[p]{r,t) + v,c[p]{r,t) (74) 

where Vext{r, t) = v{r) + u{r, t) comprises the static and td external potentials of the system 
under consideration. Apparently, Eq. ()74j) serves as a definition of an xc potential-functional 
Vxc[p]if,t) by subtracting two known potentials from the unknown KS potential-functional 
w[p]{r,t) (see Refs. js^, 38, Ig^). At least formally, everything looks as one would expect. 
But is this really the solution? Once more, the rKS formulation allows for a closer inspection 
of what we have got, because here the potential w[p]{r,t) can be given in an explicit form. 

Indeed, as the analysis of Sec. 5.B has shown, for a density p{r,t) there exists a single- 
particle td SE (Eq.EII), 

z^V(r', t) = {-IV^ + w[p] (r, t)}^(r, t) (75) 

with the td local potential (Eq. [TOj) . 

w[p]{r,t) = ^-l{Vky-k (76) 

which is determined by the density up to a purely td function {a{t)). Let us note that now 
the indefinite td function is no longer relevant, because it does not affect the resulting density. 
But it seems that Eqs. ()75I76|) lack any predictive power. They hold for any density, and one 
may wonder how the time-development of the exact density of the interacting A^-electron 
system, po{r,t), could be determined unless po is already known and used to construct 
w[po]{r,t). 

At this point it is instructive to inspect the more transparent case of static DPT. Let 
us assume for a moment that there is no second Hohenberg-Kohn theorem (HKII) and, 
thus, no variational principle. As above, however, one has a shortcut to KS-type equations 
(now applying the HKI theorem to the non-interacting system). In the rKS variant, the 
corresponding single KS-type equation can explicitely be constructed (by inserting the ansatz 
in the one-particle Schrodinger equation): 

{-|V2 + «;[p](r)}0(r) = 60(r) (77) 

Here the potential 

^[p]W = ^ + c (78) 
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is determined by the density p(r) = N(j){rY (up to a constant c). Eqs. ()77l78p show that 
any (reasonable) density is non-interacting (one-electron) v-representable. But, clearly, the 
potential-functional w[p\{r) of Eq. (|78|) as such is of no avail for determining the exact 
ground-state density, po- According to the successive steps, 

p(r) — > w[p\{r) — > (Eq.EZI) — > 0(r) — > p{r) 

any density p(r) will only reproduce itself. Obviously, the potential-functional (j78p is trivial, 
i.e. without physical meaning. In the variationally derived KS equation ()49|1 . by contrast, 
the potential-functional v^^f[p\{r) = v{r) + J[p](r) + Vxclplir) according to Eq. is of 
completely different type. The density will change in the course of the iterative solution 
of the KS equation and will (eventually) converge to the exact (or approximate) ground- 
state density Pq. Only for po, the KS equation (with v'^-l'^[pQ]{r)) will reproduce the initial 
density po. At this point, the non-trivial and the trivial potential become identical (up to a 
constant), f '^■^■^ [po] (''") = 'w[Po]{f) + c, as can be seen by comparing Eqs. (jK^ and (ffH|) . 

Yet this is not the end of the story. While hardly discussed in the literature, the triviality 
of the KS potential-functionals arising in the HKI/RGI shortcut is well-known among the 
DFT theoreticians. Rather than "defining" an xc potential-functional by subtracting the 
given one-particle potential, v{r), and the Hartree potential, J[p](r), from the trivial KS 



potential-functional, w\p 
following partioning jSy, 



r) , a non-trivial xc potential-functional can be established by the 
w[p]ir) = v^Mir) + J[p](r) + v^ir) (79) 



Here fext[p](^) is the potential-functional established by the HKI theorem for the interacting 
iV-electron system, yielding 

v{r)=v,,t[po]{r) (80) 

for the exact ground-state density, po. Letting aside the f-representability problem in the 
Vext[p] potential-functional, one can easily see that the xc potential- functional defined by 
Eq. (|79|) is the same as that arising in the variational derivation. Having established the 
existence of a non-trivial xc potential-functional, one may now use it in the KS equation in 
the familiar way, that is, together with the given one-particle potential v{r) of the considered 
system. What one gets is fully equivalent to the variationally derived result (Eqs. I49I5U|) : a 
KS equation that does not reproduce any density except for the exact ground-state density 
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Po, thereby offering the possibihty of determining po by a self-consistency procedure. The 
only difference between the variational and non-variational derivation is that in the former 
case the iterative procedure represents a well-defined search for a minimum on an energy 
surface, whereas in the non-variational approach the final step amounts to an ad hoc ansatz 
of a regula falsi type, for which the question of convergence remains open. In practice, of 
course, this distinction does not matter, because, wittingly or unwittingly, the variational 
KS derivation will serve as a safeguard for the convergence of the self-consistency cycle in 
the non-variational approach. 

Now we may come back to the time-dependent case. Lacking variationally derived KS 
equations, one can, nevertheless, establish a "non-trivial" xc potential-functional in analogy 
to Eq. (HHD: 

«;[p(t)](r-,t) =t;e.4p(t)](r,t) + J[p(t)](r,t)+t;,,[p(t)](r,t) (81) 

where Vf.xt[p{t)\{'r .,1) denotes the potential-functional established via the RGI theorem (ana- 
logue to HKI). For the exact density trajectory, po{r, t), of the interacting A^-electron system 
with the external potential Vext{fit) the RGI potential-functional gives 

Vext[po{t)]{r,i) =Vext{.r,i) (82) 

Again, we may use the non-trivial xc potential-functional w^c[p(^)](*", ^) in the (radical) KS 
equation together with Vext{f^'t)i giving rise to the td one-orbital Schrodinger equation 

i^7/^(r,t) = {-iv2 + t;,,,(r,t) + J[p(t)](r,t) + t;,,[p(t)](r,t)}^(r,t) (83) 

which is "correct" only for the exact density Po(t', t)- In the latter case the effective potential 
in Eq. (jHS)) becomes w[po](r,t), that is, the trivial KS potential-functional of Eq. (fTHjl taken 
at the exact density po(r,t). This means we have returned to our starting point and the 
original question: Will the KS Eq. (|HHj) with the potential w[pQ\{r,t) allow us to determine 
the time evolution of the exact density pQ{r,t) without knowing it beforehand? 

Let us first note that there is no longer a self-consistency cycle providing for any feedback 
to the exact solution. Eq. (jHH|) would have to be solved by time-propagation all the way 
"along" the exact density trajectory po{t). In practice, that would mean to propagate an 
"incorrect" equation, using a guessed xc potential, and starting from an approximate ground 
state density at, say, t = 0. But apart from such practical reservations, the basic question 
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is whether the propagation would succeed given the exact xc potential-functional and the 
exact density at t = 0. But even in this ideal case the answer is seen to be negative. The 
problem here is that the mapping p{t) — » w[p(t)] (Eq. 17^ is "non-local" in time (or non- 
instantaneous). What does this mean? An instantaneous potential-functional, for example, 
is the Hartree potential, J[p(t)]{r,t): the density at given time t determines the Hartree 
potential at the same moment t. But the situation is not as simple in the case of the KS 
potential-functional. This can be seen by inspecting again the KS potential-functional ()76|) 
of the rKS formulation. The first contribution on the r.h.s. of Eq. ()76p is instantaneous; 
the second term, depending on the gradient of the phase function k{r,t), requires the first 
time-derivative of the density according to the discussion in Sec. V.B; and the third term, 
being the time- derivative of the phase function k{r, t), can only be determined if the second 
time-derivative of the density is available. A similar temporal non-locality must be expected 
for Vxc[p(t)] and Vext[p(t)]- The consequence for the time- propagation of Eq. (jH^ is obvious: 
the second time-derivative of the density is not determined by the development through a 
given time t ( "past" ) so that the potential at time t is undefined unless one takes into account 
also the density trajectory beyond the point t ("future"). The inevitable conclusion is that 
the non- variational td KS equations are not suitable for predicting the time-development of 
the exact density. They would allow one to reproduce the time-development of the density 
at the orbital level, provided the td density is already given, e.g., from a solution of the full 
iV-electron td Schrodinger equation, which is of course without practical use. At a purely 
formal level of the theory, this crucial difference between predicting and merely reproducing 
is well concealed and, thus, easy to overlook. 

In view of this finding, the "causality problem" in the TDDFT linear response equations 
is to be seen from a new perspective. Rather than being the consequence of improperly 
defined xc functionals, as van Leeuwen and others have supposed, the problem seems to 
reflect the basic inadequacy of a linear response treatment for the unphysical, merely formal 
td KS equations. 

D. Linear response in the adiabatic approximation 

So far we have used the rKS concept as a tool to analyze some basic aspects of the 
TDDFT approach. Let us finally take a view at the structure of the results to be expected 



26 



at the linear response (LR) level of the theory. 

Using the adiabatic approximation for the td xc potential, both in the usual and radical 
KS versions, leads to the RPA-like equations given by Eq. (|21|). where the blocks of the 
secular matrix, A and B, are constant (cu-independent) matrices. The excitation energies, 
Um = Em — Eq, are obtained as the eigenvalues of the pseudo-eigenvalue problem 



The transition moment associated with the — m transition derives from the corresponding 
(specifically normalized) pseudo-eigenvector components according to 



The manifold of excitations obtained from these equations is determined by the configuration 
space of the secular matrix block A (note that the RPA pseudo-eigenvalues occur in pairs 
having positive and negative values, respectively). For the ordinary KS approach this means 
that the excitation manifold is that of the p-h or single excitations (with respect to the GS 
KS determinant). Here each spatial p-h configuration gives rise to 4 (primitive) spin states, 
from which one singlet and three (degenerate) triplet states can be formed. Let Uo = N/2 
and denote the number of occupied and virtual spatial KS orbital. Then the KS LR 
excitation manifold comprises ^riony solutions. The full excitation manifold of interacting 
electrons is, of course, much larger, because double and higher excitations come into play. 
It is thought that the restriction to single excitations is a consequence of the adiabatic 
approximation and the missing double and higher excitations would be accounted for by 
going beyond that approximation. Supposing that the exact td xc potential, axc[p]{r,t), 
exists, one would arrive at the same type of equations as in Eq. fl21|) . but now with u- 
dependent matrices, A{u) and B{u). In principle, this could lead to an enhanced excitation 
manifold. 

Let us now inspect the excitation manifold in the rKS case. The LR equations within 
the adiabatic approximation have the same structure (Eq. ITIj) as those of the usual KS 
approach, but there is only one occupied spatial KS orbital {no = 1). As a consequence, 
the excitation manifold comprises only excitations, that is, the excitations out of a single 
(average) KS orbital. Moreover, any spin degrees-of-freedom are missing, and, even if one 




(84) 




(85) 



a,k 
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assigns the riy spatial excitations to singlets, any triplet excitations are absent. We have 
argued that the td treatment in the rKS framework is as legitimate or not legitimate as the 
usual KS approach. This means that here the hypothetical non-adiabatic td xc potential, 
axc[p]{r , t) , must not only account for the double and higher excitations but has to restore 
already the single excitation manifold. A generation of triplet excitations appears to be 
completely impossible, because the density function p{r,t) and, thus, axclplif , t) , does not 
bear information on the spin degrees-of-freedom. Let us note that the absence of triplets 
in the rKS version does not constitute an inconsistency with the A'"-particle KS case. In 
principle, the standard KS approach too does not allow for triplets, because the perturbing 
td potential has to be local, and a local potential cannot excite triplets from a singlet 
ground state. This fact is often repressed because, in a technical sense, the usual TDDFT 
LR equations do yield triplet excitations (albeit with vanishing intensities). In the rKS 
scheme triplet excitations are neither accessible basically nor technically. 

As was noted in Sec. 4, the rKS variant is only the hmiting case of more general M- 
electron KS schemes, where the number M of non-interacting electrons may even exceed N. 
In the latter case the adiabatic approximation would produce more single excitations than 
the original interacting A^-electron system, which would mean that the non-adiabatic td xc 
potentials must eliminate spurious solutions introduced at the adiabatic level of theory. 

This shows that the well-known excitation manifold or "counting" problem of the LR 
form of TDDFT is further aggravated in the rKS (and M-electron KS) variants. Whereas it 
cannot be excluded that a hypothetical energy dependent xc potential beyond the adiabatic 
approximation might restore the single-excitation manifold of the ordinary KS scheme and, 
moreover, generate double and higher excitations, it appears more convincing to see this 
problem as an indication of the invalidity of the TDDFT equations. 

VI. SUMMARY AND CONCLUSIONS 

The title of this paper poses the question whether TDDFT is formally exact. What 
have we learned to answer that question? Let us summarize the three main topics of our 
investigation. 

First, we have observed that an error is introduced both in the td and static KS linear 
response theory if the perturbing (external) potential is given by a non-local operator. This 
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error, resulting from the neglect of the h-h and p-p matrix elements of the perturbing operator 
in the KS response equations, is of second order in the electronic repulsion. Yet for a local 
(multiplicative) potential no logical inconsistency arises, because the absence of h-h and p-p 
matrix elements does not imply a loss of information of the local operator. As stated in 
Theorem 1 of Sec. 3, this is a remarkable general property of local operators. It is a well- 
known consistency requirement of DFT that the external potentials must be local. One can 
also extend the HK and KS approach to general non-local potentials, but that would require 
the incorporation of the non-local potentials already in the definition of the HK functional, 
Fhk{p]-i and, accordingly, in the KS xc potential, Va;c[p]- As a consequence, the functionals 
would no longer be universal but depend on the respective non-local potentials. 

The problem of the non-local operators reminds us that the KS LR equations cannot be 
viewed as having an unconditional physical meaning. In fact, their validity derives from 
the underlying theory. The CPKS equations, for instance, are founded on a valid first- 
order perturbation theory for the KS equations for the perturbed A'"-electron system. In the 
td case the validity of the response equations would presuppose that the time-dependent 
extension of the KS equations be correct, that is, the td KS equations establish a formally 
exact approach to the time-development of the iV-electron density. 

Secondly, we have discussed an utmost simplification of the KS concept, referred to as 
radical Kohn-Sham (rKS) approach. Here the ground-state density is obtained from a single 
one-particle KS equation supposing a correspondingly modified xc potential. In principle, 
the rKS form of the theory is as legitimate as the usual N-particle KS approach. Whereas the 
ordinary KS approach will certainly be better suited for developing practical computational 
schemes, it is not inconceivable that the rKS variant will have some computational potential 
as well. More importantly though, the rKS formulation represents a useful pedagogical and 
analytical tool, and as such it has been used here to elucidate basic aspects of DFT and 
TDDFT. 

In the td extension of the rKS approach, having the same validity status as the ordinary td 
KS theory, a single td one-particle KS equation, though involving a possibly very complicated 
td xc potential, would allow us to determine exactly the time-development of the density 
function of the full interacting A'^-electron system, thus bypassing the A^-electron time- 
dependent Schrodinger equation. Given the richness of the phase relations of the full N- 
electron wave function, the spin degrees-of-freedom, even the permutation symmetry, all the 



29 



wealth of information seemingly absent in a one-particle orbital or the one-particle density 
function, the possibility of predicting the exact time- development at a fictitious one-particle 
level must appear fantastic. But this is what TDDFT implies. Is that expectation justified? 
Guided by the rKS approach, in the third step, we have critically reexamined the RG 
foundation of TDDFT. What has become apparent here is a phase problem corrupting the 
definitions of the RG action integral functionals, as already recognized and analyzed by van 
Leeuwen in a different way js^. Our findings fully confirm van Leeuwen's conclusion that 
the RG foundation of TDDFT, based on analogues to the HKI, HKII, and KS theorems, is 
invalid. 

But there is an alternative way of establishing the KS equations of TDDFT without 
invoking a variational principle, as formulated by Gross, van Leeuwen and others. Again, 
this idea can be analyzed at the rKS level. The mere existence of a one-particle KS type 
equation reproducing the exact td density is almost a triviality within the rKS framework. 
But due to temporal non-locality the corresponding KS potential-functional does not allow 
one to predict the time-development of the exact density of the interacting iV-electron system 
under consideration (without introducing the information of the exact density at some point). 
In its present shape, TDDFT is just a formal theory without any predictive power. As we 
have argued in Sec. 5.C, van Leeuwen's construction of Keldysh functionals does not cure 
that basic deficiency. 

Given the original RG foundation of TDDFT invalid and the design of a KS theory 
without a variational principle an illusion, it appears hardly possible to escape the conclusion 
that the idea of TDDFT, that is, the idea of a formally exact method for predicting the time- 
development of an interacting A^-electron system at the orbital level, must be abandoned. 
The TDDFT approach (in linear response form) was first introduced 25 years ago as an 
analogue to TDHF (or RPA) before any attempts at a rigorous foundation had been made. 
Without the RG theorems or another viable justification, the theory would be set back to 
the status it had in its beginning: an empirically "corrected" version of the RPA While 
TDDFT (LR) may afford an improvement over the RPA description, it cannot escape the 
RPA limitation of being an approximate method for singly excited states. 
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APPENDIX: MATRIX REPRESENTATION OF LOCAL OPERATORS 



For matrix representations of local one-particle operators the following uniqueness theorem 
holds: 

Theorem: If two local operators, v{r) and w{r), have the same p-h matrix elements with 
respect to an arbitrary partitioning of a (complete) one-particle basis into ho\e{h) and 
particle(p) states, they can differ only by a constant A, that is, w{r) = v{r) + A. 
Proof: Let (f)s{f), s = 1,2, ... , denote the functions (orbitals) of a one-particle basis and 
assume a partitioning of the orbitals such that 0i(r), . . . , 0„(r) are referred to as occupied 
or hole states and 0„+i(r), 0n+2(^)5 • • • as unoccupied or particle states. Consider two local 
operators v{r) and w{r) having the same p-h matrix elements, 

{(pa\v\(l)k) = {(f)a\w\(f)k) , k <n, a> n (A.l) 

This means that the p-h matrix elements of the difference operator, A(r) = w{r) — v{r) 
vanish: 

Aafc = {4>a\w - v\(f)k) = 0, k <n, a> n (A. 2) 

Now consider the n functions X{r)(f)i{r), I < n. These functions may be expanded in terms 
of the basis functions, 

oo 
n 

= J]AH</'fc(r), / < n (A.4) 

k=l 

yielding finite linear combinations as a consequence of Eq. ()A.2|) . Obviously, the latter 
equations can be brought to diagonal form by a suitable unitary transformation: 

X(r)^k{r) = \kk4>k{r), k <n (A.5) 

This means that all transformed diagonal matrix elements must be equal, Xkk = X, k < n, 
and the difference potential is constant: A(r) = A. 

The proof given here shows that the theorem can also stated as follows: Any local operator 
with vanishing p-h matrix elements with respect to a complete basis set and an arbitrary 
partitioning into p and h states must be a constant. An apparent objection is: What about 
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a diagonal representation of the operator? The answer is that local operators cannot be 
diagonalized properly, that is, in the Hilbert space of P functions. 

An interesting question arising in this context is if it is possible to reconstruct a local 

operator v{r) (up to a constant), if only its p-h matrix elements are given. It seems that this 
is not possible except for the special case n = 1 (one occupied state). Let us first inspect 
the case n — 1 and let 4>i{r) be the single h orbital. Expanding v{r)4>i{r) yields 

oo 

v(r)0i(r) = vuMf) + Yl (^-6) 

a=2 

where Vpq = {(f)p\v\(f)q) denote the matrix elements of v{r). Dividing this expression by (f>i{r) 
yields an explicit representation 

oo 

v(r) ^vn + Yl ^aiMr)Mr)~^ (A.7) 

a=2 

which reconstructs v{r) in terms of the p-h matrix elements Vai up to a constant, being here 
the (single) h-h matrix element, I'll. 

In obvious generalization of the case n = 1 one may proceed as follows. Let there be n 
occupied orbitals, 0i(r'), . . . ,(j)n{r),n > 1. Expanding the products, w(r)0j(r), gives 

n oo 

v{r)(f)i{r) = 'Yvki(l)k{'^)+ XI ^o-M'^)^ (A.8) 

k=l a=n+l 

where the summation on the rhs has been split into h and p parts. As in the proof above, the 
h-h block of the v matrix can be diagonalized by a suitable unitary transformation, yielding 

oo 

v{r)^i{r) = Vii^i{r) + ^ w„i0„(r), i < n (A.9) 

a=n+l 

Here Vai denote the transformed p-h matrix elements of v{r). Dividing these equations by 
the respective transformed occupied orbital, (j)i{r), leads to n different representations of 

v{r), 

oo 

v{r) ^Vii+ Va4a{r)4>i{r)~^, i < n (A.IO) 

a=n+l 

However, this does not solve the problem because the transformed p-h matrix elements, Vai 
cannot be determined without diagonalization of the h-h block of the v matrix, that is, 
without the knowledge of the h-h matrix elements. Thus, it appears that one encounters 
a similar situation as in the theoretical foundation of DFT, where the Hohenberg-Kohn 
theorem states the existence of an universial xc functional without any constructive means. 
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